A novel numerical model for protein sequences analysis based on spherical coordinates and multiple physicochemical properties of amino acids.

How to characterize short protein sequences to make an effective connection to their functions is an unsolved problem. Here we propose to map the physicochemical properties of each amino acid onto unit spheres so that each protein sequence can be represented quantitatively. We demonstrate the usefulness of this representation by applying it to the prediction of cell penetrating peptides. We show that its combination with traditional composition features yields the best performance across different datasets, among several methods compared. For the convenience of users, a web server has been established for automatic calculations of the proposed features at http://biophy.dzu.edu.cn/SNumD/.

Characteristic dysbiosis of gut microbiota of Chinese patients with diarrhea-predominant irritable bowel syndrome by an insight into the pan-microbiome.

BACKGROUND: Irritable bowel syndrome (IBS) is reported associated with the alteration of gut microbial composition termed as dysbiosis. However, the pathogenic mechanism of IBS remains unclear, while the studies of Chinese individuals are scarce. This study aimed to understand the concept of dysbiosis among patients with Chinese diarrhea-predominant IBS (IBS-D), as a degree of variance between the gut microbiomes of IBS-D population and that of a healthy population. METHODS: The patients with IBS-D were recruited (assessed according to the Rome III criteria, by IBS symptom severity score) from the Outpatient Department of Gastroenterology of Peking University Third Hospital, and volunteers as healthy controls (HCs) were enrolled, during 2013. The 16S rRNA sequences were extracted from fecal samples. Ribosomal database project resources, basic local alignment search tool, and SparCC software were used to obtain the phylotype composition of samples and the internal interactions of the microbial community. Herein, the non-parametric test, Wilcoxon rank-sum test was carried out to find the statistical significance between HC and IBS-D groups. All the P values were adjusted to q values to decrease the error rate. RESULTS: The study characterized the gut microbiomes of Chinese patients with IBS-D, and demonstrated that the dysbiosis could be characterized as directed alteration of the microbiome composition leading to greater disparity between relative abundance of two phyla, Bacteroidetes (Z = 4.77, q = 1.59 × 10) and Firmicutes (Z = -3.87, q = 5.83 × 10). Moreover, it indicated that the IBS symptom features were associated with the dysbiosis of whole gut microbiome, instead of one or several certain genera even they were dominating. Two genera, Bacteroides and Lachnospiracea incertae sedis, were identified as the core genera, meanwhile, the non-core genera contribute to a larger pan-microbiome of the gut microbiome. Furthermore, the dysbiosis in patients with IBS-D was associated with a reduction of network complexity of the interacted microbial community (HC vs. IBS-D: 639 vs. 154). The disordered metabolic functions of patients with IBS-D were identified as the potential influence of gut microbiome on the host (significant difference with q < 0.01 between HC and IBS-D). CONCLUSIONS: This study supported the view of the potential influence of gut microbiome on the symptom of Chinese patients with IBS-D, and further characterized dysbiosis in Chinese patients with IBS-D, thus provided more pathological evidences for IBS-D with the further understanding of dysbiosis.

Relationship between Gut Microbiota and Phosphorus Metabolism in Hemodialysis Patients: A Preliminary Exploration.

BACKGROUND: Hyperphosphatemia is a risk factor associated with mortality in patients on maintenance hemodialysis. Gut absorption of phosphate is the major source. Recent studies indicated that the intestinal flora of uremic patients changed a lot compared with the healthy population, and phosphorus is an essential element of bacterial survival and reproduction. The purpose of this study was to explore the role of intestinal microbiota in phosphorus metabolism. METHODS: A prospective self-control study was performed from October 2015 to January 2016. Microbial DNA was isolated from the stools of 20 healthy controls and 21 maintenance hemodialysis patients. Fourteen out of the 21 patients were treated with lanthanum carbonate for 12 weeks. Thus, stools were also collected before and after the treatment. The bacterial composition was analyzed based on 16S ribosomal RNA pyrosequencing. Bioinformatics tools, including sequence alignment, abundance profiling, and taxonomic diversity, were used in microbiome data analyses. Correlations between genera and the serum phosphorus were detected with Pearson's correlation. For visualization of the internal interactions and further measurement of the microbial community, SparCC was used to calculate the Spearman correlation coefficient with the corresponding P value between each two genera. RESULTS: Thirteen genera closely correlated with serum phosphorus and the correlation coefficient was above 0.4 (P < 0.05). We also found that 58 bacterial operational taxonomic units (OTUs) were significantly different and more decreased OTUs were identified and seven genera (P < 0.05) were obviously reduced after using the phosphate binder. Meanwhile, the microbial richness and diversity presented downward trend in hemodialysis patients compared with healthy controls and more downward trend after phosphorus reduction. The co-occurrence network of genera revealed that the network complexity of hemodialysis patients was significantly higher than that of controls, whereas treatment with lanthanum carbonate reduced the network complexity. CONCLUSIONS: Gut flora related to phosphorus metabolism in hemodialysis patients, and improving intestinal microbiota may regulate the absorption of phosphate in the intestine. The use of phosphate binder lanthanum carbonate leads to a tendency of decreasing microbial diversity and lower network complexity.

Prediction of translation initiation site for microbial genomes with TriTISA.

UNLABELLED: We report a new and simple method, TriTISA, for accurate prediction of translation initiation site (TIS) of microbial genomes. TriTISA classifies all candidate TISs into three categories based on evolutionary properties, and characterizes them in terms of Markov models. Then, it employs a Bayesian methodology for the selection of true TIS with a non-supervised, iterative procedure. Assessment on experimentally verified TIS data shows that TriTISA is overall better than all other methods of the state-of-the-art for microbial genome TIS prediction. In particular, TriTISA is shown to have a robust accuracy independent of the quality of initial annotation. AVAILABILITY: The C++ source code is freely available under the GNU GPL license via http://mech.ctb.pku.edu.cn/protisa/TriTISA.

Tree and rate estimation by local evaluation of heterochronous nucleotide data.

MOTIVATION: Heterochronous gene sequence data is important for characterizing the evolutionary processes of fast-evolving organisms such as RNA viruses. A limited set of algorithms exists for estimating the rate of nucleotide substitution and inferring phylogenetic trees from such data. The authors here present a new method, Tree and Rate Estimation by Local Evaluation (TREBLE) that robustly calculates the rate of nucleotide substitution and phylogeny with several orders of magnitude improvement in computational time. METHODS: For the basis of its rate estimation TREBLE novelly utilizes a geometric interpretation of the molecular clock assumption to deduce a local estimate of the rate of nucleotide substitution for triplets of dated sequences. Averaging the triplet estimates via a variance weighting yields a global estimate of the rate. From this value, an iterative refinement procedure relying on statistical properties of the triplets then generates a final estimate of the global rate of nucleotide substitution. The estimated global rate is then utilized to find the tree from the pairwise distance matrix via an UPGMA-like algorithm. RESULTS: Simulation studies show that TREBLE estimates the rate of nucleotide substitution with point estimates comparable with the best of available methods. Confidence intervals are comparable with that of BEAST. TREBLE's phylogenetic reconstruction is significantly improved over the other distance matrix method but not as accurate as the Bayesian algorithm. Compared with three other algorithms, TREBLE reduces computational time by a minimum factor of 3000. Relative to the algorithm with the most accurate estimates for the rate of nucleotide substitution (i.e. BEAST), TREBLE is over 10,000 times more computationally efficient. AVAILABILITY: jdobrien.bol.ucla.edu/TREBLE.html

Accuracy improvement for identifying translation initiation sites in microbial genomes.

MOTIVATION: At present the computational gene identification methods in microbial genomes have a high prediction accuracy of verified translation termination site (3' end), but a much lower accuracy of the translation initiation site (TIS, 5' end). The latter is important to the analysis and the understanding of the putative protein of a gene and the regulatory machinery of the translation. Improving the accuracy of prediction of TIS is one of the remaining open problems. RESULTS: In this paper, we develop a four-component statistical model to describe the TIS of prokaryotic genes. The model incorporates several features with biological meanings, including the correlation between translation termination site and TIS of genes, the sequence content around the start codon; the sequence content of the consensus signal related to ribosomal binding sites (RBSs), and the correlation between TIS and the upstream consensus signal. An entirely non-supervised training system is constructed, which takes as input a set of annotated coding open reading frames (ORFs) by any gene finder, and gives as output a set of organism-specific parameters (without any prior knowledge or empirical constants and formulas). The novel algorithm is tested on a set of reliable datasets of genes from Escherichia coli and Bacillus subtillis. MED-Start may correctly predict 95.4% of the start sites of 195 experimentally confirmed E.coli genes, 96.6% of 58 reliable B.subtillis genes. Moreover, the test results indicate that the algorithm gives higher accuracy for more reliable datasets, and is robust to the variation of gene length. MED-Start may be used as a postprocessor for a gene finder. After processing by our program, the improvement of gene start prediction of gene finder system is remarkable, e.g. the accuracy of TIS predicted by MED 1.0 increases from 61.7 to 91.5% for 854 E.coli verified genes, while that by GLIMMER 2.02 increases from 63.2 to 92.0% for the same dataset. These results show that our algorithm is one of the most accurate methods to identify TIS of prokaryotic genomes. AVAILABILITY: The program MED-Start can be accessed through the website of CTB at Peking University: http://ctb.pku.edu.cn/main/SheGroup/MED_Start.htm.